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ABSTRACT 


Intensities of auroral and magnetospher ic hiss generated by the 
Cherenkov radiation process of electrons in the lower magnetosphere 
are calculated with respect to a realistic model of the earth* s magneto- 
sphere. In this calculation, the magnetic field is expressed by the 
'*Mead-Fairfield Model" (EOS, 53^, 1099, 1972), and a static model of the 
iono-magnetospheric plasma distribution is constructed by accumulated 
data obtained by recent satellite observations (Alouette-I, -II, ISIS-I, 
0G0-4, -6 and Explorer 22). The energy range of hiss producing electrons 
and the frequency range of produced VLF in the computation are 100 eV 
to 200 keV, and 2 to 200 kHz, respectively. Intensities with a maxi- 
mum around 10 kHz, of the order of 10 watts m”^ Hz observed at 
satellite altitudes and of 10" ^ watts m~^ Hz ^ at the ground seem to 
be ascribale to the incoherent Cherenkov emission from soft electrons 
with a differential energy spectrum E"^ having an intensity of the order 
of 10^ cm"^ sec"-*- str - ^ eV"-*- at 100 eV. It is shown that the frequency 

of hiss spectral density peak at geomagnetic latitude 70° night side 

o 

is around 20 kHz for the soft spectrum (~ E ) electrons, which shifts 
to 10 kHz for the hard spectrum (~ E“ • ) electrons. The maximum hiss 
intensity produced by soft electrons is more than one order higher than 
that of hard electron produced hiss. Higher rate of hiss occurrence 
in the daytime side, particularly in the soft electron precipitation zone 
in the morning sector, and less association of auroral hiss in nighttime 
sectors must be, therefore, due to the local time dependence of the 
energy spectra of precipitating electrons rather than the difference in 
the geomagnetic field and in the geoplasma distributions. 


// 



1 . INTRODUCTION 


Estimations of hiss intensity associated with auroras have been made 
recently by Lim and Laaspere (1972) , in which the contribution of the 
relatively low energy electrons below 1 KeV is taken into consideration* 

Their calculation has shown the importance of low energy electrons on 
auroral hiss generation which has been disregarded in previous theoretical 
works, mainly due to the lack of observational information of low energy 
electrons and of reliability of theoretical formula (McKenzie, 1963, 1967; 
Mansfield, 1964, 1967; Liemohn, 1965; Jorgensen, 1968). Measurements of 
low energy auroral electrons by means of rockets (Evans, 1968; Ogilvie , 

1968; Feldman et al., 1971) and of satellites (Frank and Ackerson, 1972, 
Hoffman, 1969; Burch, 1968) have provided abundant information. Comparisons of 
low energy electron spectra observed by satellites and rockets as well 
as their variations with latitude and time have been shown recently in 
great detail, indicating the dynamics of these electrons in the earth 1 s 
magnetosphere (Rees, 1969; Hones et al. 1971). Furthermore, the better 
correlation of auroral hiss occurrence with precipitations of low energy 
e lectrons rather than those of high energy electrons has been shown by 
the latest analyses of simultaneous satellite observations of auroral 
VLF-hiss and of electron precipitation events (Hoffman and Laaspere, 

1972, Gurnett and Frank, 1972). 

On the other hand, the calculation by Lim and Laaspere (1972) is 
limited for the daytime polar region (78° geomagnetic latitude) where 
the plasma distribution is derived from the theoretical model of Banks 
and Holzer (1969) and the magnetic field is assumed to be that of the 
centered dipole. Obviously the earth 1 s magnetic field differs significantly 
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from that of the centered dipole particularly near the so-called cusp 
region where most auroral daytime events including VLF-hiss take place. 

Owing to numerous satellite observations, the distribution of plasma, 

1. e., densities of thermal electrons, protons and ions, in the earth's 
ionosphere and magnetosphere has been known in fair detail. Since 
these plasma data are essential for the estimation of hiss-generation 
and the propagation through these media a more realistic model of 

the earth's plasma distribution than that used by Lim-Laaspere (1972) 
is worth constructing, with indications of its variation with latitude, 
altitude and local time. 

The purpose of this paper is, therefore, to recalculate the auroral 
hiss intensity and intensities of VLF-hiss in the magnetosphere produced 
from magnetospheric low energy electrons, by making use of a realistic 
model of the earth's magnetic field and of plasma distributions. The 
latter is constructed by making use of the accumulated satellite data 
of electron, proton and ion density distribution^ obtained by Alouette I, 

II, ISIS-I, 0G0-4 and 5. The geomagnetic field is expressed by the 
"Mead-Fairf ield 1972 model " ( EOS , 53, 1097, 1972). 

2. FORMULATIONS OF THE CALCULATION 

A charged particle moving along its helical trajectory in the tenuous 
magnetoplasma emits Cherenkov radiation in addition to the synchrotron rad- 
iation, if its velocity parallel to the static magnetic field is larger 
than the local phase velocity of the wave in the plasma. In contrast to 
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synchrotron radiation which is emitted as a result of decelerating gyro- 
motion of the particle due to the static magnetic field imbedded in the 
plasma, Cherenkov radiation is emitted from the polarized wake produced 
in the medium a fast moving electric charge. 

The frequency range of the wave associated with these polarized 
wakes of plasma medium is wide but limited below a certain frequency 
characterized by the plasma and magnetic field, as shown in the next 
section. Propagations of these waves in the earth's magnetosphere are 
very slow, dispersive and mostly along the field line until their fre- 
quencies approach another characteristic frequency of the magnetoplasma, 
called the lower hybrid resonance (LHR) frequency, f^. Since the propa- 
gation of the wave in the magnetoplasma above the LHR is essentially due 
to the mobility of electrons in the medium, which is far larger than 
that of ions, the polarization of the wave is the same sense as that of 
electron-gyration, i . e ., right-handed with respect to the direction of 
the magnetic field. This mode of wave propagation is called the whistler 
mode (R-mode VLF) . 

The phase velocity of waves propagating in a magnetoplasma is shown 
in Figure 1, as a function of wave angular frequency a)(=2nf) for the 
case when the electron gyrofrequency f^ and the plasma frequency f^ are 
10 kHz and 30 kHz, respectively. As can be seen from, this figure, 
phase velocities . of waves in the whistler mode are far slower than .the 
light velocity c, and decrease with wave frequency, approaching that 
of Alfven waves below the ion- gyrofrequency , Q^, of the plasma. 
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It should be noted that the upper limit of whistler mode frequency 
is fjj if f H < f n , in contrast to the case shown in Figure 1. 

2.1. The Formula of Cherenkov Radiation in Magnetoplasma 

The Cherenkov emission from an electron moving along its helical 
trajectory in the magnetic field imbedded in a cold, collisionless plasma 
has been described in many theoretical papers. Although the final ex- 
pressions are more or less the same, the derivation of the formula of 
Cherenkov radiation from an incident spiraling electron in plasma is 
not unique . 

Eidman (1958) has applied the so-called Hamiltonian method known 
in the field theory of quantum electrodynamics (Heither, 1944) for the 
estimate of emission rate from the polarized anisotropic plasma medium. 
This is further extended to the relativistic form by Liemohn (1965). 

McKenzie (1963, 1967) and Mansfield (1964, 1967) have used classic 
forms of the solution of Maxwell *s equations for the polarized electric 
field and the current in the anisotropic plasma induced by a spiraling 
incident electron by means of Fourier transforms. According to Trulsen 
and Fejer (1970), Eidman 1 s results are more accurate in general than 
Liemohn* s, even though his derivation is based on the non-dispersive 
medium with non-re lativistic formulation. 

The present calculation follows essentially the Mansfield method 
as has been used by Jorgensen (1968) and by Lim-Laaspere (1972). Since 
Mansfield*s notations and definitions of plasma parameters are different 
from those commonly used in the physics of plasmas and of the ionosphere. 
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the basic formulae are rewritten in the Appendix B: notations, definitions 

and constants are listed in Appendix A, which are mostly taken from 
books by Stix (1962) and by Ratcliffe (1962). 

As shown in Appendix B, by (B-50) , the intensity of Cherenkov 
emission per unit time, i.e., the power of Cherenkov radiation, P G (in 
watt/cm^) from an electron with kinetic energy, E (in keV) averaged 
over the spiraling period Tj|(=2rr/r^) is expressed as functions of the 
radiated wave frequency f (=uu/2rr) , of the magnetic field intensity, B 
(in gauss) and of plasma parameters. The plasma parameters are repre- 
sented by a plasma frequency f^(=2n7n o ) and the magnetic field is 
expressed by the gyrofrequency of electron, f^ (=Qe/2TT). 

The dependence of P Q on and is expressed by the inverse di- 
electric tensor (cosQ) shown by (B-39), where the analytic expressions 
of S, P and D as functions of t h and of can be seen from (B-13). 

Most VLF measurements are made not on the total intensity over the 
whole observable frequency range but rather on the differential intensity, 
i.e., on intensities at different frequencies.. Therefore, the power per 
frequency, dP Q /df (in watts/Hz el.) is a more useful presentation for the 
comparison with experimental results. This can be written from (B-52) 
as follows : 



( 2 . 1 ) 
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where J Q , J-^ are the Bessel functions of the zero-th and the first order, 
respectively with the argument 
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j 


“> 01 . . 

v sm© 

Qe J 


( 2 . 2 ) 


The so-called Cherenkov angle 0 Q , in (2.2) is given by the following 
relativistic formula ; 


where 


and 


COS0 O - — 

HI J 

j3j_ - j3 coso?, /3 1 | = fi sina 

Q m /e(E +2 trigC 2 ) 

E + tn e c 2 


(2.3) 


(2.4) 


E and n^c^ are the kinetic energy (in KeV) , and the rest mass of electron 

o 

m^c = 510.98 (in KeV), respectively, a is the pitch angle of helical 
trajectory of the electron with respect to the static magnetic field 
direction. 

The two values of the refractive index of a plasma medium stand for 


n 


j 


2 


-b„ -(-l)j A 0 

2S 


(2.5) 


with 


A 0 


b„ = 


p_: s 
18 1 


+ D 2 - S 2 - PS 


c D = ^ — 9 (S 2 - D 2 - PS) + P (S z - D z ) 


?2 _ t,2 ' 


18 1 


( 2 . 6 ) 


and 


2 


( 2 . 8 ) 



2.2. The Dispersion Relation of Waves in the Plasma 

As shown in Figure 1, the propagation velocity of waves in the 
magnetoplasma changes with the angle (0) of propagation vector |k of 
the wave from the direction of magnetic field tB. The directional 
dependence of propagation velocity further changes with wave frequency. 
These angular and frequency dependences of wave propagation in a plasma 
can be found by calculating the dispersion relation of the wave in the 
plasma. 

Since the Cherenkov radiation occurs only when the velocity of an 
impinging electron along the field line, V ^ , is larger than the phase 
velocity, Vph> °f the wave in the plasma, the refractive index n (^c/vph) 
is one of the essential parameters in the present calculation. The ratio 
of the critical velocity of the electron, V c (below which no Cherenkov 
emission takes place in the plasma) to the velocity of light in vacuum, 
j3 c = V c /c, is therefore equal to the inverse of the refractive index 
of the medium. This can be seen from the Cherenkov condition (B-28) 
with 9=0. The /3-value of the electron is shown in Figure 2 as a function 
of the kinetic energy of the electron E. The right side scale of 

ordinates indicates the value of the refractive index n corresponding to 
the critical value of /3 C and the kinetic energy E of the electron. It 
should be noted that in the energy range below around 100 keV, jS is 
proportional to the square-root of E as can be expressed by the non- 
relativistic expression of kinetic energy, i.e. 
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On the other hand, in the energy range above roughly 2 MeV, the 
velocity of the electron can be regarded as almost the velocity of light 
in vacuum, where the approximation called "the extreme relativistic, 

~ l* 1 ' can be used. It is in this energy domain where the Cherenkov 
radiation has been used for the detection of charged particle in high 
energy physics (Oostens et. al, 1967; Yuang et al., 1969). 

In plasma, however, the propagation of the polarization in the 
medium, which is the electromagnetic wave propagation in the medium, 
is so slow that the Cherenkov emission is important rather in the non- 
relativistic energies than in the extreme relativistic energy domains. 

As shown in the following sections, the rate of Cherenkov emission in 
the plasma is larger at lower energies than at higher energies, as long 
as the velocity of incident electrons is larger than that of plasma elec- 
trons. However, if the energy of the incident electron approaches the 
thermal energy of the medium, the formula of Cherenkov emission based on 
the cold plasma is no longer correct. The correction is, however, very 
small as long as the energy of the incident electron is above the thermal 
energy of the plasma as shown recently by Trulsen and Fejer (1970) and 
by Trulsen (1971). 

The dispersion relation of cold plasma is derived from the de- 
terminant of the dielectric tensor of the plasma, i.e., 

|t(cos9) | = 0 

where TC (cos9) is given by (B-16)*. As can be seen from (B-17) with 
(B-18), this leads to the fourth order equation of n. Therefore, there 

o 

are two solutions for the real values of n . Using the matrix elements 
shown in (B-16), the solution of (B-17) can be written in the following 
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form 

9 B + F 
n z = — = — 
2A 


(2.9) 


where 

F 2 = (RL - PS) sin 4 © + 4P 2 D 2 cos 2 © (2.10) 

A, B, S, P,, D, R and L are given by (B-18). 

Equation (2.9) presents one of the dispersion relations of the plasma 
and is shown in Figure 3, by plotting n^ against ou/Qe* 

Substitution of (B-18) into (2.9), brings the following form of 
the dispersion relation 


2 q = (.. n2 -L) 

(Sn 2 -RL) (n 2 -P) 


( 2 . 11 ) 


From (2.11) we find the following propagating wave modes in the 
plasma ; 

(i) Propagation along the magnetic field, by setting 9=0, we get 
P=0 , n 2 = R and n 2 = L (2.12) 


^Footnote: Since the propagation vector k is taken in ZY-plane, the 

matrix elements in T (cos9) are different from those of Stix 
(Eq. 20 in p. 11, 1962), in which k is taken in ZX-plane 
while the direction B is the same in Z-axis. Because of 
the axial symmetry of the plasma characteristics around Z- 
axis , the dispersion relation (2.9) is the same as that of 
Stix (Eq. 21 in p. 11, 1962). 
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(ii) Propagation across the magnetic field, by setting 9=7, 

n^ = and n^ = P (2.13) 

S, P, R and L are analytic functions of u), fij, Ilj and 9 as shown 
in (B-13) , where suffix j stands for electron (j=l), proton (j=2) and 
other ions in the plasma (j = 3, 4 ). 

In the actual earth’s ionosphere and magnetosphere, the gyrofrequency 
fie is not necessarily smaller than the plasma frequency TI e as shown in 
Figure 1. In Figure 3, therefore, the dispersion relations are shown 
for three cases i.e. (a) < n e , (b) fi e = IT e and (c) fie > n e , 

assuming fi e = 836 kHz* One of the common features in these figures 
which will be called the "dispersion diagram" of the plasma, is that wave 
modes in the diagrams consist of two basic domains, distinguished by the 
refraction index (n) of the plasma, i.e., the one (A) with n^ ^ 1 
and the other (B) with n^ > 1 . 

In the domain (A) , there are three cutoff frequencies which are 
given by setting n^ = 0 in (B-17) • S^nce C=RPL, these cutoff frequencies 
correspond to R=0, P=0 and L=0 modes, respectively. The mode R=0 , which 
is right-handed in circular polarization, corresponds to X-mode in Figure 1. 
P=0 and L=0 modes are left-handed in circular polarization and correspond 
to the 0- and Z-modes, respectively. The nomenclature X- , 0- and Z-modes 
are made for the classification of radio waves in the ionospht_*^ 5 based 
on the so-called Appleton-Hartree dispersion formula (Ratcliffe, ?962). 

The Appleton-Hartrie formula can be obtained from Equation (2.9) by 


making an approximation 
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Qj « cju < Q e 

In the domain (B), there are at least four resonant frequencies in 
which the index of refraction diverges, i.e., n^ co } and the phase 
velocity of the wave approaches to zero, Vp^ =0. As shown in Appendix 
C, these resonant frequencies are the upper hybrid frequency, the 

electron gyrofrequency , the plasma frequency n e and the lower hybrid 
frequency For the case of only an electron plasma, as used for drawing 

Figure 3, cju^h ”* 0. In the multi -component plasma, Qj and Ilj appear also as 
resonance ^frequencies . In Figure 3, the right-handed and left-handed 
modes are indicated by dotted and hatched domains respectively. 

As can be seen from these figures (Fig. 1, 2 and 3) the waves produced 
by the Cherenkov process are limited in this domain, (i.e., n^ > 1 or 
v ph < c) . The propagation modes and the directional phase velocities of 
the wave, particularly those around the resonance frequencies are dis- 
cussed in detail in Appendix C. 


2.3 The Plasma Data 


Since the rotation of the electric field vector of a right-handed cir- 
cular polarized wave (R-mode wave) is the same as the gyration of an electron 
in a magnetic field, the waves produced by the motion of electrons in the 
magnetoplasma are R-mode waves (in brief, R-waves) . 

As shown in Figure 1, the frequency range of R-waves is limited be- 


tween two resonance frequencies, and u) LH - Due to other resonance 

frequencies, Qe a nd II e , which exist between uo^ and uo LK , these waves are 


divided into two modes, i.e., the one is called the whistler mode and 
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the other is the R-mode EM wave in a narrow frequency band. If 
as shown in Figure 1 and Figure 3(a), the whistler mode is limited be- 
tween Q e (where the propagation is parallel to B, i.e., 9=0) and 
(where the propagation is perpendicular to B, i.e., 9=tt/2). The other 
R-mode wave frequency uo is limited in n e - ti) - w h ere the propagation 

is parallel to iB at cd = II e and perpendicular to iB (i.e., 9 = ^*) at u) = 

The changes of propagation direction with wave frequency at a given point 
in the magnetosphere are shown by calculating the refractive index surface 
as a. function of frequency, f (for example, Jorgensen, 1968, Fig. 9). 

If Qe > n e , then as shown in Figure 3(c), the frequency range of the 
whistler mode is oUyj ^ uu ^ II e , where the propagation is longitudinal* 
(9—0) at cjo — n e and transversal* (9 — tt/ 2) near oj s cjo^p: as before. 

The other R-mode wave is in a narrow band of II e ^ U) ^ (jutjh with trans- 
verse propagation near and longitudinal near n e . The band shrinks 
further into the narrow band Qe < (ju < U)^ at larger values of the re- 
fractive index, n ^ > 50 or for slower phase velocities, Vp^ < 10 km/sec. 
(Vph < 0.3c). This narrow band of R-mode is connected to the whistler 
mode, if n e = Qe where the propagation is longitudinal (9—0 propagation) 
in both modes. This case is shown in Figure 3(b). 

It should be noted that the most part of auroral hiss is in the 
whistler mode with a) < (if Q e < TI e ) or oj < TI e (if II e < Qe). 

* Footnote: Here the words longitudinal and transverse are used for 

parallel and perpendicular to the imposed magnetic field, 
as commonly used in the ionosphere physics (Ratcliffe, 1962) 
In plasma physics, however, they are used for jk ||E 0 and 
ki£ 0 while parallel, perpendicular are used for k ||B 0 and 

k jlBq, respectively (Stix, 1962, p.39). 
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Due to the increasing magnetic field intensity B with decreasing 
altitudes, the whistler mode VLF produced at high altitudes in the earth* s 
magnetosphere approaches the local lower hybrid frequency uu^ as it 
propagates toward the ionosphere. Correspondingly, the direction of 
propagation becomes more transversal as the wave approaches the iono- 
sphere from its original longitudinal propagation at large distances. 
Consequently, some of the lower frequency VLF -waves are reflected 
back into the magnetosphere where the propagation becomes more longi- 
tudinal again until the wave enters into the other hemisphere. Repeating 
these process, VLF -waves of whistler mode can be trapped between two 
hemispheres and in the earth's magnetosphere. An important role of these 
trapped VLF-waves in reducing the lifetime of electrons in the radiation 
belt, i.e., the formation of the "slot" around L — 3, has been discussed 
recently by Russel and Thorne (1970). 

In this respect, the calculation of flg , Tl 0 , uu^ and as functions 

of the altitude, of the latitude and of the local time are rather im- 
portant not only for the calculation of VLF-production and its propagation 
in the magnetosphere but also for the study of dynamics of trapped charged 
particles in the earth's radiation belts. 

2.3 The Model of the Geomagnetic Field 

As discussed in the introduction, the earth's magnetic field has 
been assumed to be that of the centered dipole in the previous theoretical 
works on VLF emission and its propagation in the magnetosphere (Helliwell, 
1965; Jorgensen, 1968; Lim-Laaspere , 1972). The deviation of the actual 
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field configuration from the centered dipole field can be seen from 
Figure 4, in which several magnetic field lines are drawn for comparison. 
The dipole field is shown by thin lines and the field calculated by a 
model of the earth's magnetic field called "MF- 72 (Mead -Fairfield 1972 
model) " is indicated by heavy lines. The field lines starting from the 
daytime side of the earth's surface and from the nighttime side are dis- 
tinguished by full lines and dashed lines, respectively. The MF-72 is 
constructed by least squares fitting to magnetic field measurements 
obtained by latest available satellites. Explorer 33, 34, 41 and 43 
(Mead and Fairfield, 1972). This model has been elaborated by separating 
the data and the corresponding parameters in the model into two cases of 
geomagnetic conditions, i.e., for quiet times (Kp < 2) and the other for 
disturbed periods (Kp ^ 2). Since the auroral hiss occurs during 
polar and magne tospheric substorms, MF-72 model for Kp > 2 conditions 
should be the logical choice for the purpose of the present calculation. 

The difference between the two cases is, however, a small affect on 
the result of the present calculation, although the difference of these 
models from the centered dipole field is significant as can be seen from 
the comparison shown in Figure 4. The difference between the two cases 
(quiet and disturbed) of MF-72 is significant- at large distances from the 
earth (larger than 10 R^) while most of the VLF-hiss production takes place 
within the domain of distance less than 10 Rg . Furthermore, the magneto- 
spheric model should depend on the solar wind and on the interplanetary 
magnetic field (IMF). The effects of these variations are rather complex 
functions of time and space to be considered in the present subject. This 
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requires a wide range of parameters hard to fit into a single model. In 
Figure 5, the field lines and corresponding iso-gyrof requency lines of 
electrons in the Kp < 2 model used in the present calculation are 
shown, where (a) is for geocentric distances less than 15 Rg and (b) for 
distances less than 5 Rg , respectively. The units between the marks 
shown along field lines are 1.0 Rg and 0.2 Rg , respectively. The unit of 
gyrof requency is kHz. These figures correspond to Helliwell*s dipole 
field charts (Helliwell, 1965, pp. 316-319). 

2.5. The Distributions of Magnetospheric Plasma 

As discussed Appendix C, the upper and the lower hybrid resonance 
frequencies, and are functions of gyro- and plasma - frequencies 

of charged particles in the earth* s magnetoplasma. In order to compute 
these characteristic frequencies, a model of the spatial distributions 1 
of electrons, protons and 0 + - ions in the magnetosphere (including the 
ionosphere) has been built by using the latest available data obtained 
by the following satellites and investigators (indicated in the parenthesis); 

Alouette-I and II (Chan and Colin, 1969; Norton, 1969; Jackson, 1967, 

1969; Rycroft and Thomas, 1970), 0G0-4 and 6 (Mayr et al., 19 70; 

Bnnton et al., 1971; Taylor, 1972; Taylor and Walsh, 1972) and Explorer 
22 (Brace et al., 1970). 

Since these data do not cover whole domains of the magnetosphere, 
extrapolations as well as interpolations of the data have been made in 
several parts in the model, by utilizing the theoretical works on the 
plasma transport in the polar-magnetosphere (Banks and Holzer, 1969 a,b; 

Banks et al . , 1971). Some parts of the results are shown in Figures 6, 
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and 7 by plotting the number densities of electrons, protons and 0 + 
ions against geomagnetic latitudes for the daytime and nighttime sides, re- 
spectively*. Since the extremely disturbed periods are not included in 
these data, this geoplasma distribution can be regarded as a model for 
Kp < 2. 

2.6. Characteristic Frequencies in the Magnetosphere. 

Based on the model of the MF-72 geomagnetic field and that of the 
iono-magnetospheric plasma distributions described in the previous 
sections, the upper and the lower hybrid resonance frequencies, f^ 
and f^ are calculated as functions of the altitude, latitude and local 
time. These are shown in Figures 8 and 9, together with the electron 
gyrofrequency , f^, and the plasma frequency, f^, plotting them against 
the distance from the ground measured along the magnetic field lines 
(in km) for different geomagnetic latitudes; (a) 60°, (b) 70° and (c) 80° 
for daytime, and for (a) 60° > (b) 70°, (c) 75° and (d) 80° for night 
time, respectively, where the electron gyrofrequency = Q e /2n, the 

^Footnote: Most of original data below 1000 km level are given as a 

function of geographic latitude, while most magnetospheric data are 
given against dip or geomagnetic latitudes. Since the difference, 
particularly the longitudinal dependence due to the difference of 
these coordinate system is not critical to the present calculations 
in which only midday and midnight conditions are considered, all 
parameters are shown against geomagnetic latitudes. The difference 
between dip and geomagnetic latitude is important only near the 
polar lower ionospheric regions. 
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plasma frequency f N = n 0 /2TT, the upper hybrid resonance frequency 
fyn = CJL) uh^ tt anc * t * ie l° wer hybrid resonance frequency f LH = cjo^r/ 217 are 
calculated by the equations (A-l) , (A-3) , (C-15) and (C-16), respectively. 
Some significant aspects of these figures are; (i) there is an alti- 
tude range at each latitude both in the day side and in the night side, 
where the plasma frequency is smaller than the electron gyrofrequency , 
i.e., f N < f H , corresponding to the dispersion diagram Figure 3(c), while 
outside of this altitude range, the dispersion relation corresponds to 
the diagrams of Figure 3(a) and Figure 1 (f^ > f^.) . (ii) There is a 

minimum of f around the latitude of 70° X for each altitude, corre- 
L>n 

sponding to the so-called ion-trough (Taylor, 1972; Taylor and Walsh, 
1972), although f T decreases with the latitude in general (Brice and 
Smith, 1965). (iii) There is a maximum of f^ along each field line on 
both the day side and night side, forming a barrier for the propagation of 
VLF of frequencies lower than 10 kHz from the outer magnetosphere to the 
lower ionosphere. (iv) In the night, due to the low electron density 
below the altitudes around F-maximum, (below 350 km level) there is 
another barrier for high frequency side of VLF above several hundreds 
kHz as can be seen from Figure 9. 

It would be worth noting that the physical meaning of f UH and f^ : 
can be described as follows; under the high frequency waves, the 

motion of positive ions is negligible as compared to that of electrons. 

The electrons, on the other hand, undergo oscillatory motion by the 
periodic force of the wave which resonate at the plasma frequency. In 
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the magnetoplasma, the motion of electrons becomes circular due to the 

static field and their gyrating motions are imposed on the plasma 

o 2 k 

oscillation. The resulting frequency is therefore ud^ = (0^ 4- n o ) 2 

which approaches the plasma frequency as the magnetic field is reduced 
to zero, i.e., -* II 0 as B -> 0. 

Contrary to the high frequency disturbance, at the very low fre- 
quencies, electrons can be regarded as a qua si -homogeneous background 
tightly bound to the static magnetic field, while oscillatory motions 
of ions under slow periodic field variations wiLl resonate at a certain 
frequency that is cr^jin a manner that electrons oscillate at 

3. CALCULATIONS AND THE RESULTS 

3.1. Energy Dependence of Hiss -Production 

With respect to the plasma distribution described in the previous 
sections, the power of Cherenkov emission per frequency (watts/Hz) by 
an electron are calculated by the equation (2.1). along the field lines 
in the magnetosphere. These powers are shown in Figures 10 and 11 for the 
daytime side and for the nighttime side, respectively. Since these 
differential powers (watts/Hz) of Cherenkov emission depend on the 
energy of the hiss-producing electrons, figures are shown for (a) 

E = 100 keV, (b) E = 1 keV for dayside (Figure 10), and (a) E - 100 keV, 
(b) E = 10 keV, (c) E = 1 keV and (d) E = 0 . 1 keV for nightside (Figure 
11), respectively. These results correspond to tho^ shown in Figure 10 
of Jorgensen’s paper (Jorgensen, 1968). The present results indicate 
the following; (i) As shown in the calculation by Lim and Laaspere (1972), 
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the power of Cherenkov emission, contributing the VLF-hiss in the magneto- 
sphere increases with decreasing energy of electrons. (ii) The power 
of emission per electron is larger for higher frequencies for a given energy 
of electron, but (iii) this increase is not monotonic, particularly when 
the total emission per frequency in the magnetosphere is considered. 

This will be discussed in the next section. 

3.2. Distributions of Hiss-Production in the Magnetosphere 

The differential intensity of hiss production by electrons per unit 
volume in the magnetosphere, DP/Df (in watts Hz cm 8 ) can be shown as a 
function of s, the arc-length from the earth's surface along the geo- 
magnetic field line (in km), by integrating (2.1) with respect to the 
pitch angle, a, of electrons and to their energy, E. Assuming the iso- 
tropic pitch angfe distribution, this is given by 

2tt N e (E) • sin a da dE (^) . (3.1) 

where dP/df (in watts/Hz .el . ) is given by (2.1), N e (E) = j e (E)/v e is 

3 -1 -1 

the number density of electrons (cm - str keV ). The 
velocity v e = /^ c (in cm sec -1 ) of the electron with kinetic 
energy, E (in keV) is given by (2.4) and is shown in Fig. 2. 

Using the energy spectrum of electrons observed by Hoffman and 
Laaspere (1972), i.e. j e (E) = 4 x 10 8 • E -2 (cm -2 sec -1 str -1 keV -1 ) , 
the integration of (3.1) with respect to the electron-energy E is per- 
formed numerically. This energy spectrum of electron has been also 
applied by Lim and Laaspere (1972) in their calculation. The results 
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for nightside are shown in Figure 12, (a) for f = 100, 10 kHz and (b) 

for f = 50, 20 and 5 kHz, respectively, where the lower limit of the energy 

integral, E c is assumed to be 0.1 keV. 

3.3. Variation of Hiss Production with Energy Spectrum of Electrons 

The total intensity of auroral hiss per frequency produced by 
electrons in the magnetosphere can be estimated by integrating (3.1) 
with respect to the effective volume of production. Assuming the per- 
fect ducting of hiss propagation between the magnetosphere and the 
ground level for the frequency range shown in Figures 8 (day side) and 
9 (night side), this is given by 

D(P/f) = J S (g|) • S(s) • ds (3.2) 

s o 

where DP/Df is given by (3.1) and shown in Figures 12 (a) and (b), 

the cross- sectional area S(s), (in m z ) at the distance s, which 

corresponds to the unit area S Q at the ground, s Q , is given by 

the magnetic flux conservation, i.e., 

“ <5(7)> s o * 0.3) 

B Q and B(s) are the total intensity of geomagnetic field at s Q 
and s, respectively, and S Q = 1 m^. 

In order to see the effect of the energy spectrum of electrons on the 
rate of hiss production, the similar calculation with a harder spectrum 
j e (E) = 8.10^ E~^*^ (cm"^ sec“^ str~^ keV ^) is made with respect to the 
same geomagnetoplasma model. The total intensity of electron is normalized 
at E c = 1. keV for both spectra, i.e., the total intensity above 1 keV 
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is the same for both cases, J(E C )=4 x 10^ (els/cm^ sec str). The re- 
sults are shown in Figure 13 (a,b) with a full line and a dashed line 
for the hard and soft spectra of electrons, where (a) stands for 80° X 
dayside and (b) for 70° X nightside respectively. Notice that the 
vertical scale for the hard spectrum case (full line) is shown in the 
right side. 

This figure indicates the following: (1) The peak intensity of 

hiss produced by electrons of soft energy spectrum is more than one 
order higher for the hiss produced by soft spectrum electrons. (2) The 
wave frequency of peak intensity is also higher for the hiss generated 
by soft spectrum electrons than the one by hard spectrum electrons. 

These results agree qualitatively with those obtained by the recent 
calculation by Lim and Laaspere (1972). However, the spectral density 
has a peak intensity around 20 kHz, in contrast to 70 kHz in Lim-Laaspere 1 s 
calculation, although the spectrum of electron energy is the same. The 
peak intensity of the present result for the soft electrons is roughly half 
that of Lim-Laaspere . This difference seems to be reasonable because of 
the shift of the peak toward lower frequency where the power of Cherenkov 
emission per electron is generally less than at higher frequency as shown 
xn Figures 10 and 11. 

As mentioned in the Intrpduction , the geomagneto -plasma model used 
by Lim and Laaspere (1972) is less accurate than the one used in the 
present calculation, especially at great distances from the earth where 
most of the low frequency hiss below 30 kHz is produced. Since in these 
calculations, the perfect ducting along the field line is assumed for 
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the propagation of VLF-hiss in the magnetosphere, the peak intensity 
corresponds roughly to the maximum volume of hiss production in the 
magnetosphere. As can be seen from Figures 8 and 9, therefore, the 
frequency of peak intensity should be around 20 kHz for both dayside 
and nightside. 

The frequency of the peak intensity of hiss produced by hard spectrum 
electrons, 10 kHz, is rather in good agreement with Jorgensen* s cal- 
culation (1968), although the energy spectrum of electrons is not exactly 
the same. The shift of the peak intensity to lower frequency with the 
harder energy spectrum of hiss producing electrons has been also stated 
in Lim-Laaspere r s paper (1972). As can be seen from Figures 10 and 11, 
the increasing rate of Vl*F-emission with decreasing energy of electrons 
is larger at higher frequency. The surpressing effect of tte harder spectrum 
on hiss production is, therefore, relatively dominant at higher frequencies. 
The shift of peak frequency of VLF hiss tcward lower frequency by the harder 
energy spectrum of source electrons can be explained by this effect. 

3.4. Discussions 

As presented in Jorgensen *s paper (1968), the peak intensity of 
auroral hiss observed at the ground (Byrd station, Antarctica) and by the 
0G0-2 satellite is both around 10 kHz (which is by no means constant but 
mostly below 30 kHz). In. this respect, the models of the geomagnetic 
field and of the geoplasma distribution used in the present calculation 
seem to be very close to reality, giving proper values of the frequency 
for the peak intensity and of spectral density of auroral hiss. 
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On the other hand, as in the previous work by Jorgensen (1968) 
and by Lim and Laaspere (1972), the perfect ducting is also assumed in 
the present calculation for the propagation of all VLF with frequencies 
between the upper limit (electron gyrofrequency or plasma-frequency) 
and the lower limit (lower hybrid frequency). The conditions and the 
effect of the ducting for the magnetosphe ric hiss propagation have been 
discussed recently by Manoranjan Rao et al. (1973), who pointed out that 
near LHR- frequencies , the propagation of hiss requires large density 
enhancements along the field line instead of the previously assumed a 
few percent deviation for the ducting. These will be investigated further 
at well as the contribution of low energy electrons below 100 eV . It 
is known that the first effect reduces the intensity near LHR and the 
latter effect increases the peak intensity and the frequency of maximum 
intensities at all latitudes. 
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APPENDIX A 


Notations and 

Definitions 

a: 

pitch angle of electrons 

jB(=¥/c) , 

where V is the velocity of electrons, km/sec. 

$-\\ > 

are normal and parallel component of J3 


with respect to the magnetic field line IB. 

B: 

magnetic induction, gauss 

e o : 

dielectric permittivity of plasma 

E: 

kinetic energy of electron, keV 

f : 

frequency of wave, sec”"** 

n ( =v ph/ c )> 

refractive index of the plasma, where Vp^ 


is the phase velocity of the wave. 

T(=l/f) 

period of the wave, sec. 

f H (=ne/2TT) , 

gyrofrequency of electron, sec~^. 

f N( =n o/ 2r 0 , 

plasma frequency, sec“^. 


~ 1,78 x 10 ^ B, the angular gyrofrequency 
of electron, sec“ 1 where B is in gauss. (A-l) 

n e (=/S^-) 

* e o “e 

~ 5.65 x 10^ /N e , electron plasma angular 
frequency, sec ^ where N e is the electron 
density, cm . (A-2) 

"o " We 2 + n p 

O O — 1 

+ 11^ ) 2 , angular plasma frequency, sec . (A-3) 

n. = / e2N i z i , 

angular plasma frequency of i-th ion, sec ^ 

V e D mi 

where and m^ are the number density. 


atomic number and the rest mass of the i-th 


ion. (A-4) 

Zj e B 

fij. = 1 

m i c 

angular gyrofrequency of the i-th ion, sec"^. 
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A-2 


s : 


S(s) : 



j e ( E > : 


d£o 

df 


DP . 
df 


D(P/f): 


length of the geomagnetic field-line from 
the earth's surface, km. 

cross section of the geomagnetic tube at a 

distance s, corresponding to the unit area 
o 

at s=o, nr. 

2 

plasma current, coulomb /cm sec 

differential energy spectrum of electrons, 

-2 -1 -1 -1 
cm sec str keV 

time-average of radiated power over a gyro- 
period of electron, watts per electron, 
radiated power averaged over a period of 
gyration, per electron per radiation-frequency, 
watts/Hz - electron. 

differential intensity of hiss radiated from 
unit volume of the earth 1 s magnetosphere 
(Cherenkov emissivity of the magnetosphere 
per radiation frequency), watts/Hz cm . 
the total intensity of hiss per unit frequency 
produced by electrons in the magnetosphere, 
watts/Hz . 


(»>) 


Physical constants 
Light velocity in vacuum. 
Rest mass of the electron, 
Rest mass of the proton, 
Electronic charge, 


c = 

2 . 997 ^ x 10 5 

km/ sec 

“e 

= 5 .110-^ x 10 2 

keV/c 2 

“p 

= 9.3826 x 10 5 

keV/c 2 

e = 

■ 1.602 1 x 10“ 19 
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APPENDIX, B 


Electric field E in. the plasma is related to the source current, j s , 
produced by incident particles through Maxwell's equation (in MKS units) 
as follows; 


and 


dD 

V X ® = dt J s 


9D = 
dt 



\V x E = 


c© 

dt 


(B-l) 

(B-2) 


(B-3) 


where e c is the dielectric permittivity of the plasma, jp is the plasma 
current, which arises from polarization of the plasma by jj g , and 
can be expressed by the summation of individual currents, jj of 
each species of the plasma, i.e.. 


where 


h = S j J 

J 


h ' “J E 


(8-4) 


<S 7 j is the conductivity tensor of the plasma which is obtained by 
solving the Lorentz force equation for individual species of the 
plasma , 


qj( E 


(r,t) 


+ Vj (r,t) x B) = mj 


dVj (r ,t) 
dt 


(B-5) 


where m. , Va and OB are the rest mass, the electric charge, the 

j J J 

/ velocity of the j-th specie and the static magnetic field in the 

plasma, respectively. 

Since the current jj can be written as 


jj = qj Nj Vj 


2G 



B-2 

where Nj is the number density of the j-th specie of the plasma, Fourier 
transform of (B-5) gives 

- i W Jj + 6j Qj X J. = e o n. 2 . E (B-6) 


where Qj and II j are the angular gyro-frequency and the angular plasma- 
frequency given by 


«j Oj = q j B 


na * 


and II . = 
J 



(B-7) 


respective ly . 

Comparison of (B-6) with jj = cfj j , we get the following express- 
ion for the conductivity tensor as a function of Q j , rr j , and uu, 



(B-8) 


Fourier transform of (B-l) , (B-2) and (B-3) yields 


'k x B = - cue 0 + i j s 

where 

K = I - -Hi ... 

iou e c 


(B-9) 


(BIO) 


with 


( sr T = e «j . 

L j J 
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B-3 


and n is the unit tensor, 


H si 



Introducing the propagation vector m = tk 


(B-9) can be written 


[(inn - n^IL) + K] • E = — i j< 


0jp o 


whe re 



(B- 11) 


(B- 12) 


is called the dielectric tensor of the plasma (Stix, 1962), and the 
elements are given as functions of wave- frequency , uu, gyro-f requency , Q j , 
and plasma- frequency II j , i.e. 


S = 1 - £ 


V 


j OJ 2 - fij2 


D = E 6 j n j‘ 


and 


j ou(tu 2 - n 2) 

p = 1 - Ej n j 2 

J CU 2 


Sj = + i (B-13) 


If the direction of the static magnetic field, B, is taken in z-axis, 
the character of the plasma has the axial symmetry around the z-axis. 
Without losing the generality, the propagation vector of a wave, /k 9 can 
be taken in yz-plane, then 


k x - 0, and ik = ayk sin 9 + a z k cos 9 
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(B-14) 



B-4 


where 0 is the angle between B and k , 


a x , a y , a z and k x , ky , k z are unit vectors and wave numbers in 
x, y, z directions, respectively. 

Substitution of (B-12) into (B-10) yields 


TEE = 


1 


a) e, 


- J i 


(B-15) 


where 


S-n^ 


- iD 


TC(cos 0) =( iD S-n 2 C os 2 0 n2 sin Q cos e 


? 

rr- sin 


0 cos 0 P-n 2 sin 2 0 


(B-16) 


and 


k - 


a n -1 


B +C 


(B-17) 


where 


A = S sin 2 0 + P cos 2 0 
B = PS(1 + cos 2 ©) + RL sin 2 0 
C = PRL, R = S + D and L = S - D 


(B-18) 


1=0 gives the dispersion relation of waves in the collisionless, 
cold plasma. 

From (B-15) , we get 

= T _1 (cos 0) J s Ck,co) (B-19) 
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B-5 

By the Fourier integral theorem (Fourier's inversion formula for (B-19)), 
E(r,tj = f f £(k,u ;) • e 1 ^ “ k * r > dk dcju 

we get 

E(|r,t) = |- JJ T' 1 (cos 0) J s (k,uj) • e 1 ^ UJt k dk ^ (B-20) 

O -QO 

The source current j s can be expressed by 

<j s (r) = q s -\V s (t) • 6 (r - r s (t» (B-21) 

The Fourier transform of this yields 


j s (|k,(ju) = 


<2tt)‘ 


J # ”\V s (t) • exp [ i (ik - or s - out)] dt (B-22) 


Since a particle spirals around z-axis, moving along (B direction with 
the velocity, Vs > the radial vector of the particle, or s (t), is given by 


r s (t) = a x r Q cos •Q / ' t + a r Q sin C2 / -t + a z V||t 


(B-23) 


and 


\V s (t) = -a x V x sin t + a V x cos Q / t + a z V| 


(B-24) 


where r G = — i- is the gyro-radius of spiraling particle 

Cl' = ^s/y and y is the Lorentz factor of the particle, i.e.. 
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B-6 


Since k = ,n V = j3c , the product of (B-14) and (B-23) can be 

written as 

ik-r s - cut = § sin Q 7 t+ (nj3||cos0 - 1) cut 


where 



sinQ 


The integrand of (B-22) is therefore written in the following form 


j s (k,w) = -^TT f (-a x V x sin O' t + a y V ± cos Cl' t + a z V||t) 

• exp [i(5 sin 0 7 t + (n/3 cosQ - 1) cut] • dt 

(B - 2 5 ) 

n 00 rt oo S 

- T^TTT 2 J (ia x V x J s (?) + a V S C§) + a 2 V„J s (?)) 

(2tt) s= -00 -00 ^ 

x exp [i(s Q 7 + cun jSjj cos0 - uj)t] • dt 


The second expression of (B-25) is derived by using the following 
formula of Bessel function; 


and 


sin ^ s ^ n ^ 

cos Q'fe 1 ? sin n ' t 


00 , 

i ^ J s (5) 

S= -00 
°° S 

s r J s «) 

S= -co b 


e is n't 
e is n't 


i I sin n't 


^ J s (§) 

S= -oo 


e is n't 
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B - 7 


Using the Fourier integral representation of delta function, i.e.. 


«<•> 'sf' 1 


cut 


dt 


(B-25) can be written in the following form 


T (Jc,cu) = Js s (ia x V x -J s '(S) + a V x 7 J s (§) + a z V„J s (§)) 
b (2rr)° s= -co ^ * 


x 6 ( s Q 7 + cun j3|| cos0 - u)) 


(B-26) 


This expression indicates that the radiation occurs only when the 
source current is in the following condition, 


s Q 7 + uo n /3|| cosQ - oo = 0 


(B-27) 


In the case of s=0, this gives the Cherenkov condition. 


cos© = — i— 

0n n 


(B-28) 


The total instantaneous radiation from a spiraling particle is given 


by 


P (ir s it) 


= (F • d ^s 


- q s (E + W s x B)* v s 
= q s E(ir s ,t). V s (t) 


(B-29) 
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B-8 


Since 


d 


Ik = k dk dfi 


t x?- up 

( TT~ 


■) 


dn dQ, 


dQ = dcp*d (cosG) 


and 


5(nu>/3|| cosO + sfif-ou) 


1 

ujn^|| 


6 (cos© 


gj - sftj 
uan /3|[ 


By substituting (B-26) into (B-20) , (B-29) can be written in the following 
expression. 


P(r s > ^ ) 


iq 


e o(2rr) 


f ° n dn f 

*-oo 


\ dou • Q s («r s ,t) dQ 
c J 


(B-30) 


where 


1 00 / 
Qs (^s.t) =1 (cosG) £ V (?) • 6(cos9 - — ^ sQ ) 

— - b ojnpji 


S= -oo 


-1 § sin Cl't + i(l-nj8|| cos©) cot , (B-31) 

V s (c s ,t) 


and 


'V s (§) = ia x V ± J s '(I) + a yVi f J s (§) + a z V|,J s (|) (B . 32) 

Due to the 6- function in the integrand, the result of integration of (B-30) 
by dQ can be obtained by replacing § and cos© by and cos0 g , respectively, 
where 





n* sin9 g 


and cos9 s 



n &w 


(B-33) 
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B-9 

The instantaneous radiation power from a spiraling charged particle 
is therefore given by. 


iq 4 00 CO co 

P(rs,t) = So"(2ttc)3 L n dn J ^ 2tt 2 Q s / (° : s >n 2 ,t) (B-34) 

' V —00 S = “OO 


where 


Q S / ( ir s» n2 » t ) ='T" 1 (cos9 s ) S V s *(?) * e -:L(m " s) nt • V s (?) (B-35) 

m= -o° 


The time-dependence of (B-34) can be eliminated by taking the time- 
average of P(tr s ,t) over a revolution-period of the particle, 2tt/c/ , 


POr s> = §7 J^? r s»t) dt 


As can be seen by (B-35), the integration of P( o: s ,t) by t over a period 
2tt/c/ vanishes except for m = s, due to the term exp [-i(m-s) fl / t]. Thus 
we get 


iq 2 


P(rs) = e“(2nc)3 J^° dn P*^ 277 .? _ QT^sT^T (B-36) 


00 S= -00 


where 


Qs / < r s» n 2 ) = V s *(?) • T-^cosQg) - ¥ s (§) 
\V S *(? S ) = (-ia x V x J s '(§), a y V x | J s ( 0 , a z V||J s (§>) 

34 


(B-37) 


(B-38) 



B-10 


TC~ 1 (cos9) = 


T(cos0) 


11 

it: 12 

-iti 3 

21 

t 22 

c 23 

31 

t 32 

t 33 


and 


v s (5) - 


/ia x VjLJs # (5) 



(B-39) 


(B-40) 


Elements of 1 ■*'(cos9 s ) in (B-39), t^j's are calculated from the inverse 
matrix of (B-16). The result is as follows: 


t n - SP - An 2 , A = S sin 2 9 + P cos 2 9 

t 12 = " n 2 sin 2 9) 

t ]_3 = Dn 2 sin9 • cos9 

t 2 1 = *‘12 

t 2 2 = n 4 sin 2 9 - n 2 (P + S • sin 2 9) + SP (B-41) 

t 23 = ( n2 " s ) • n 2 • sin9 • cos9 

fc 31 = _t 13 
fc 32 = fc 23> a nd 

fc 33 = c os 2 9 - n 2 S* (1 4- cos 2 9) + (D 2 - S 2 ) 
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B-ll 


From (B-16) , the determinant ofT(cos9 s ) can be written into the 
following form. 


|'IC(cos9 g ) | = S*n^ + b s *n^ + c s 

= S' (a 2 - n^) • (n 2 - n 2 2 ) 


where 


(B-42) 


n -^2 


~ b s + A s 
2 • S 



2 • S 


(B-43) 


A s = (b s 2 - 4'S.c s )^ 

b s = * C ■ 5) tt') 2 - (P - S) + D 2 - s 2 - PS 


(B-44) 


(B-45) 


and 


c« - 


^n 2 


(1 - ou Q^) 2 


(S 2 - 


D Z - PS) + P • (S' 


- D 2 ) 


(B-46) 


These expressions indicate that the integrand of (B-36) is regular 

except at four roots of (B-42), n = + n-^ and + n. 2 . The integration of 
o 

(B-36) by n can be therefore done by means of the residue theorem with 
respect to two poles of n^, i.e., n^=n^ j , (j=l and 2) # With manipulations 
of three matrixes products, n -integration of (B-36) yields, 


p(r s ) 


2 2 00 
S 3 J | uo | • duj 

j=l 4TTe 0 c 2 V|| 


£ Qs ' ( c s > n j 2 ) 


S= -oo 


(B-47) 
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B-12 


where 


Q s '(r s , nj 2 )- [tuV^CJs ') 2 + t^V ^) 2 Jg 2 + t 33 V ll 2 J s 2 
+ ti2^J. 2 ' “ * ' 2J S J s 

+ t 23 V x V|| • J s 2 ]/ S'(-l)j(n 2 2 - n-^ 2 ) 


As can be seen from (B-33) , the argument of Bessel function depends on 
two parameters s and j, i.e.. 


§ 





sin9 s 


where the angle of wave propagation with respect to a static magnetic 
field, 0 g , is given instead of (B-27) by discrete values of s and j, 

CM9 .,J - a - 5 >/nj 0, (B 

The result for s=0 in the above expressions, corresponds to Cherenkov 
radiation as shown by (B-28). The total power of Cherenkov radiation 
from a spiraling particle in the plasma is, therefore, given by setting 
s=0 in the above expressions. . - - - * 
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B-13 


2n , 2, 


P 0 = E fa, da, J Q 2 (Sj) + ^ * fc ll < n j > (5j> 

° U j — l “°° J 

+ 2/3 J .^t 13 (n,2)J 0 (5 j )J 1 « j )]/S(-l) 3 (>>2 2 - »1 2 ) (B - 50) 


where Jq (^j) 's are replaced by using the relation 

Jq(x)= -Jj^x), and 

U)0x 


h = 


Q' J 


tis Sin9_, cos9„ = 


011 n i 


(B-51) 


The power per frequency, d? 0 y/d£, is more convenient for comparisons 
with the experimental results than the total power of emission over whole 
frequency range. Thus, from (B-51), we obtain the following expression 
of Cherenkov radiation per emission frequency, 




fc 33 (*j 2 ) [J 0 2 (lj) + (|i) 2 J l 2 (5j) 


t ll( n j 2 > 
fc 33 < n j 2 ) 


2 ( 2*0 

Pll 


«j> J 1 




fc 13 < n j 2 ) 
c 33 ( n j 2 ) 


(B-52) 


where 


Aq =( b o 2 - 4 Sc o> % 

b = P ' 1 + D 2 - s 2 - PS 
° 


(B-53) 


(B-54) 
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B-14 


c 


o 


S 2 -P 2 -PS 

j3|i 2 


+ P(S 2 - D 2 ) 


(B-55) 


t u ( n . 2 ) = -rij 2 . S + SP + ~~ 

p ir 

t 33<"/) * - S) - JL + (D 2 . S 2j 


(B-56) 


fc 13( n j 2 ) = J^Z JafP II 2 - 1 


In these expressions, S, P, and D are the function of angular frequency 
of the wave in a plasma, which are given by (B-41) . 
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APPENDIX C 


First, let us consider the case of a two component plasma which con- 
sists of equal numbers of electrons and protons. 

Putting j=e for electrons and j=p for protons in (B-13) and 9=0 in 
(B-18), we get the following expression which corresponds to ( 2 . 12 ), 


and 


where 


n ,, 2 = ^ - <*±-0* On - n J 
(uj - Qg) (cju + Qp) 

n|| 2 = ci ) 2 + a>Q e - Og Qp - n 0 2 

(uu + n e ) (a; - Qp) 


n2 = j2 + n2 

1A o e T u p 


, for n =R (C-l) 


, for n^-L (C - 2 ) 


(C-3) 


In these derivations, terms of Qp are neglected as compared to those of 

o , ^ 2 
or and Q 0 . 

Solutions of numerators of (C-l), and (C-2), give the cutoff 
frequencies 


and 


n * 

» R ~ Pe + =2- 

u e 


n z 

OUp ~ o 


, for R-mode (C-4) 


fL 


, for L-mode (C-5) 


For the P-mode, P=0 gives 0 ^,= H e . (C- 6 ) 

Equations (C-l) and (C-2) indicate also that uu=C^ and au=Qp are resonant 
frequencies of' the wave propagating parallel to the magnetic field, i.e. 


2 

nil -* °° 
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For the wave propagating perpendicular to the field, the resonance 

TT 

frequencies are obtained from (2*13), with 9= ^ 

",2 - ^ - “ r 2 ^ 2 - » L 2 ) Eor „ 2= (c . 

(a; - (o>2 - 


n x 2 = 


(JU 2 - n 2 
O 


where cu^ and u) L are cutoff frequencies of R and L-modes, given by (C-4) 
and (C-5) , respectively. 

(C-7) indicates o)^ and uu^ are the resonance frequencies of the 
wave propagating perpendicular to the magnetic field and called the upper 
and the lower hybrid frequencies, respectively. 

As can be seen from Eq. (2*9), these are given by A=0 , which is 
written further by (B-18) 

tg 2 6 = - | (C-9) 


For 9- 2 ", this requires S=0, from which we get 


" 4 ' + “p 2 + n o 2 > + [“eV + n a 2(n p 2 + £yy] . 0 

The solution gives, 


(ju 2 = °e 2 + Q p 2 + n o 2 li + I - 4(Q e 2 n p 2 '+ n e 2 0p 2 + n^Oeflp) 

2 / (Qe 2 + ^p 2 + n o 2 )? 

^e^ n o \ + fl ~ 2 ^e^^p^ ^e^p 

2 \ (^e 2 + n o 2 ) 2 



3 


Taking the plus sign in the bracket, we get 


to. 


UH 


2 ss + n o 2 


(C-10) 


and, for the minus sign in the bracket, we get 

uj lh 2 - Q e 2 ^p 2 + n e 2 Hp 2 + n e 2 o e n p 
+ n D 2 ~~ 


By making use of = ^p^e » this can be written as, 


“LH l fi p 2 + n p2 + 




(C-ll) 


In the case of the multi-component plasma, S=0 gives the following formula, 


IL 


hn 2 + hj 2 


+ “P + “J + ... - 1 

w 2 - Pg 2 to 2 - n^2 u)2 _ Q.2 


(C-12) 


Since a2 . ^,2 fe' 


where 


e = 1 


5.10 


-4 


l p 1836 


(C-12) can be approximated by 


where 



uj 4 - uu 2 (n 2 + Q e 2 ) + n e 2 (n 2 


n 2 = n e 2 + n p 2 + n 2 + 


*) = i 

- n e 2 ) = o 


(C-13) 


(C-14) 
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Similarly to (C-10) and (C-ll), from (C-13), we get 


U) 2 ^ q 2 . yr 2 

UH e 11 


(C-15) 


and 


By defining 


^LH = 


2 ^ Q e 2 (n 2 - n e 2 ) 


n 2 + Q e 2 


1 = 2 Q?i m e z i 


(C - 1 6 ) 


(C-17) 




^eff i 

where is the fraction of i-th ion species in the total position ion. 
(016) can be approximated by (Brice and Smith 1964) 


JL 

WLH 2 


= Meff 




(C-18) 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2 . 


Figure 3. 


The phase velocity, Vp^ in km/ sec of waves in plasma as a 
function of angular frequency of the wave, cju in sec"'*’ in 
various propagation modes in the plasma, where the 
gyrofrequency of electron, f^ and the electron-plasma 
frequency f N are assumed to be 10 kHz and 30 kHz, respectively. 
Other notations are described in the list of Appendix A. 

The critical velocity of electron /J=V/c as a function of 
kinetic energy E (in KeV) , and the refractive index n=c/Vph, 
of medium, corresponding to the critical velocity. The value 
of n is indicated in the right side of ordinates. The top 
scale of abscissa is the kinetic energy in eV for ; given 
Square of refractive index n^ of the plasma given by (2.9) 
versus uo/C^ where kHz is assumed for computation. The 

dotted and hatched domains indicate possible propagation modes 
of waves for the right-handed and left-handed polarization, 
respectively. The numbers attached to boundaries of each 
domain are angle of propagation direction from the direction 
of magnetic field in degree. Dashed lines in each domains 
indicate 45° propagation from the field line. u>jjh stands for 
the upper hybrid frequency given by (C-10). (a), (b) and (c) 

in Figure 3 correspond for the cases of TI e =2Q e , II e =Q e and 
n e =0.5 Qg , respectively 
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Figure 4 


Figure 5 


Figure 6 


Figure 7 


Comparison of the Mead-Fair field model-1972 (MF-72, Kp < 2) 
earth* s magnetic field (heavy lines) with the centered dipole 
field (thin lines). The field lines are shown in the noon- 
midnight meridional plane where the field lines of MF-72 
starting from the daytime side of the surface and from the 
nighttime side are drawn by full lines and dashed lines, 
respectively. 

Field lines, arc length and electron gyrofrequency of the 
Mead-Fairf ield 1972 (MF-72, Kp < 2) model of. the earth*s 
magnetic field. The units of gyrofrequency is kHz. The in- 
terval between marks along field lines are 1.0 Rg in (a) and 
0.2 R e in (b) , respectively. 

Spatial distributions of the plasma constituents in the day- 
time side of the ionosphere and magnetosphere. The number 
density (in cm“^) of electrons (e~ ), protons (H + ) and oxygen 
ions (0+) are plotted against geomagnetic latitude by dashed 
lines, by full lines and by dash-dotted lines, respectively. 
Figures are shown for the altitudes of (a) 500 km, (b) 4000 km 
above the earth* s surface and (c) 5 Rg geocentric distance 
(1.9 x 10^ km above the surface). 

Spatial distribution of the geoplasma in nighttime side. No- 
tations of lines are the same as those in Figure 6, and the 
altitudes are (a) 300 km, (b) 600 km, (c) 1000 km, (d) 4000 km, 
(e) 3 R E and (f) 10 R E , respectively. 
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Figure 8. 


Figure 9. 


Figure 10. 


Figure 11. 


The altitude dependence of the characteristic frequencies of 
magnetoplasma in the daytime side of the earth f s ionosphere 
and magnetosphere, i.e., the upper hybrid resonance frequency 
(UHR) , the gyrofrequency of electron, f H , the plasma 

frequency, f^ and the lower hybrid resonance frequency (LHR) , 
^LH* The altitude is measured along the field line from the 
earth* s surface, s (in km) for the geomagnetic latitude at 
the surface (a) 80°, (b) 70° and (c) 60°, respectively. 

The same. as Figure 8 for the nighttime side, along the field 
line from the geomagnetic latitude (a) 80°, (b) 75°, (c) 70° 
and (d) 60°, respectively. 

The power generated at Cherenkov radiation by an electron with 
energy E (in keV) and pitch angle 0° at frequencies between 
1 kHz and 1 MHz as a function of distance from the ground along 
the field line of geomagnetic latitude 80° in daytime side. 

Two figures correspond to the energy of electron (a) 100 keV 
and (b) 1 keV. 

The same as Figure 10 except for different energies of electrons 
along the field line of 70° geomagnetic latitude from the 
ground in the nighttime side. Each figure corresponds to the 
energy of electrons, (a) 100 keV, (b) 10 keV, (c) 1. keV, and 
(d) 0.1 keV, respectively. 
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Figure 12. The power of Cherenkov emission per unit volume in the mag- 
netosphere along 70° field line by the electrons with an 
energy spectrum j e (E) = 4 x 10^ E~^ (cm“^ sec"^ str“^ keV’^). 

(a) for f=100 and 10 kHz, and (b) for f=50, 20, 5 and 2 kHz, 
respectively. 

Figure 13. The variation of Cherenkov hiss intensity with frequency of the 
emitted wave and its dependence on the energy spectrum of the 
hiss producing electrons in night side magnetosphere along 
field line (a) from 80° geomagnetic latitude dayside and (b) 
from 70° geomagnetic latitude nightside , respectively. The 
vertical scale for the soft spectrum electrons, E”^ (dashed line) 
and for the hard spectrum electrons, (full line) are 

in left side and right side, respectively. 
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